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Abstract 

We study the quark angular momentum distribution in the nucleon 
within a light-front covariant quark model. Special emphasis is put into 
the orbital angular momentum: a quantity which is very sensitive to the 
relativistic treatment of the spin in a light-front dynamical approach. 
Discrepancies with the predictions of the low-energy traditional quark 
models where relativistic spin effects are neglected, are visible also after 
perturbative evolution to higher momentum scales. Orbital angular 
momentum distributions and their contribution to the spin sum rule are 
calculated for different phenomenological mass operators and compared 
with the results of the MIT bag model. 
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1 Introduction 



The EMC measurement of the integrated helicity parton distributions [|IJ] (for 
recent results at SLAC see @) triggered the interest in a deeper understand- 
ing of how the total angular momentum of the nucleon is shared among its 
constituents. It was found that the fraction of the nucleon spin carried by 
the quarks was rather small, at variance with the most naive quark model 
expectation where the proton spin is (almost) entirely built from the spin of 
the quarks. Among the different explanations of these discrepancies, it was 
proposed || that the (overlooked until then) polarization of the gluons might 
also contribute to the singlet axial charge through the axial anomaly. In that 
case, experimental data would be compatible with a rather large fraction of the 
spin carried by the quarks (AS = 0.45 ± 0.09 in a recent world data analysis 
Q), though still far away from the non-relativistic quark model predictions. 

One of the most important issues raised by the 'spin crisis' was the need 
for considering all the possible sources of angular momentum in the nucleon. 
Therefore, the spin sum rule should read || ||: 

\mQ 2 ) + ^g{Q 2 ) + L q(Q 2 ) + L g (Q 2 ) = \, (l) 

where |AS(Q 2 ) (Ag(Q 2 )) is the spin carried by the quarks and antiquarks 
(gluons) and L q {Q 2 ) (L g (Q 2 )) is the orbital angular momentum (OAM) con- 
tribution of the quarks (gluons) J7|]. 

The significant role of OAM was pointed out several years ago |5|, but 
the problem was rigorously formulated only recently, when a gauge invariant 
definition of the quark and gluon (twist-two) operators was proposed |9|, [10]. 
Besides, there has been a big effort to derive evolution equations (at one- 
loop level) for OAM observables || [11]. At the present there is only one 
gauge invariant definition of quark OAM || with known Q 2 evolution and 
experimentally accessible (for a discussion cfr. Ref. (12|). Such a definition 



holds for reference frames with definite nucleon polarization and the OAM 
distribution could be measured through the forward limit of skewed parton 
distributions. 

One peculiar feature, already expected by general arguments |jL H], is ex- 
plicitly realized in the evolution equations, namely that the OAM distribution 
is coupled to the helicity parton distributions. As a consequence, OAM contri- 
butions can be generated, through evolution at higher scales, even in the case 
of vanishing OAM components at low initial hadronic scale. This is indeed the 
case for most hadronic models where the quarks are arranged, in the ground 
state, in a I = S-wave configuration, stressing the crucial role of Q 2 evolution 
for the evaluation of spin observables or, in other words, the roughness of the 
identification of constituent quarks with partons at all energy scales. 
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From a quantitative point of view, some studies are currently available 



fT3| , [T3| , p~5f| . In particular, in ref. ||13|| , the OAM distributions have been cal- 
culated for a number of hadronic models. As a first step these quantities are 
evaluated at the hadronic (low energy) scale and then evolved to the experi- 
mental Q 2 >> /ig scale by using the Leading Order (LO) evolution procedure 
recently established ||. One central conclusion of that work is that a sizeable 
initial OAM distributions can deeply influence the final high-energy results. 
As a consequence, a clear difference arises between non-relativistic and rela- 
tivistic models: while the former usually give a tiny OAM contribution at /Xq, 
the latter may give rise to sizeable effects at high x that persist after evolution. 
To this respect, OAM distributions are useful quantities to assess the relevance 
of relativistic effects in the hadronic models of the nucleon. 

In a recent study we investigated the consequences of a light-front treat- 
ment of relativistic spin effects on the helicity distributions [16|, [17] and in the 
present paper we want to enlarge our analysis to OAM investigating in detail 
the predictions of a light-front covariant quark model. As a matter of fact, the 
spin dynamics can be discussed within the light-front approach in a way which 
respects covariance requirements and particularly suitable to discuss deep in- 
elastic polarized processes, both at the hadronic [|15| , 16 , |T7"H and high-energy 
(partonic) scale |16|, [L7]]. We will show that light-front covariant quark mod- 
els (LFCQM) predict a non-vanishing OAM distribution whose main features 
survive after evolution. We will also see that these predictions hold for a vari- 
ety of mass operators indicating that the relevant ingredient is the relativistic 
treatment of the spin wave functions, absent in many traditional formulations 
of the quark model. The comparison with other relativistic models (MIT bag 
model) and the analysis of the moments that enter the spin sum rule will allow 
us to assess the reliability of LFCQM. 



2 OAM at the hadronic scale 

In the recent years a quark model-based approach has been developed for com- 



puting the non-perturbative inputs in the evolution equations |18| and describ- 
ing polarized and unpolarized parton distributions. Schematically, the central 
assumption is that at some low-energy scale (/ig) the nucleon is made up of 
valence quarks that can be identified with the constituents of the quark model 
(or the bag model in alternative treatments). Therefore, the non-perturbative 
boundary conditions can be evaluated by using low-energy models of the nu- 
cleon. Subsequent refinements led to the inclusion of non-perturbative gluons 
and sea at the hadronic scale /ig [l9[, as well as the explicit partonic content 
of the constituent quarks 



By following such a procedure we will assume that at the hadronic scale 
only valence quarks are resolved so that the quark helicity distribution g±(x, nl) 
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for a given flavour a is given in terms of the momentum density of the valence 
quarks: 



where x is the Bjorken variable, the mass of the nucleon and k + is defined 
as a function of the parton momentum as k + = \/k 2 + m 2 + k z . The polarized 
momentum densities are defined as: 

3 1 + rx (i) - - 
ni\k) = (N, J z = i/ 2 | £p a — ^(fc, - k)\N, J z = 1/2) , (3) 

i=l 

where V a is the flavour projector. 



An analogous definition can be worked out for the OAM distributions JL3 



where the density of the angular momentum is defined in the usual way: 



L z (k) = (N, J z = 1/2| ]T -i(ki x )-AK - k)\N, J z = 1/2) . (5) 



i=l 



From the previous definitions one can recover the result L z (x,/j,q) = 
obtained assuming S-wave quarks only, and non-relativistic approximation (i.e. 
L z (k) = 0). A more complicated example within non-relativistic dynamics, 
is given by models whose nucleon wave function is a superposition of various 



SU(6) components, such as the Isgur-Karl model [^TJ . The non-vanishing 
contribution to L z {k) in these cases is due to the D-state (or higher waves) 
admixture. For example, in the Isgur-Karl model, the OAM distribution of 
Eq. (|) results to be proportional to the D-State probability a 2 D |]T3| and its 
contribution is very small (oq = —0.067). 

This situation is radically changed in a light-front covariant quark model. 
In light-front dynamics (LFD) p2[ , the specific partition of the Poincare alge- 
bra into kinetic and Hamiltonian generators leads to several simplifications of 
the relativistic many-body problem such as, for example, the clean separation 
of the center of mass motion. The prize to pay is that the description of an- 
gular momentum is rather complicated. Not all the generators of rotations, in 
fact, belong to the kinetic subgroup, and hence the angular momentum oper- 
ator is, in general, interaction dependent. 

For this reason, in the phenomenological applications of LFD to the quark 
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model, it is customary to work in the Bakamjian-Thomas construction [p2 



that is adding a phenomenological interaction to the free mass operator, only. 
However, the resulting total angular momentum operator, although interaction 
free, does not satisfy ordinary composition rules. In order to restore them, a 
unitary transformation of the Hilbert space, known as Melosh Rotation (MR), 
has to be performed. In particular, if the nucleon is in a S-wave state, such 
rotation acts only on the spin part of the wave function. 
The D 1 / 2 representation of the MR is given by: 

, = ( m + . + M-«-(ix£,) 

{(m + uj + k z ) 2 + kl) 1 / 2 

As a result, motion and spin are now intimately correlated as it is required 
by a relativistic theory. The MR can be interpreted as the boost transforma- 
tions required to move from the rest frame of each subsystem (quark) to the 
rest frame of the total system (nucleon). 

In the present study we will not investigate SU(6) breaking effects in 
spin-isospin space: the nucleon wave function will correspond to a S-wave. 
This simplifying assumption ensures that the non-vanishing OAM contribu- 
tion originates from pure relativistic effects due to the treatment of the spin 
in light-front dynamics. Indeed the MR gives rise to a non-vanishing angular 
momentum density even if the spatial wave function corresponds to a S-wave 
and the angular momentum density can be written, for a SU (6) symmetric 
spin-isospin wave function, as 

L z (k) = \- kl n(k ) , (7) 

3 (m + uj + k z ) 2 + kl 

where n(k ) is the total momentum density, defined in the usual way : 

3 

n(k) = (8) 



i=\ 



and normalized to the number of particles (j n(k)dk = 3). Recalling the 



expressions for the polarized densities that enter the helicity distributions [16 



17 : 



nl{k)-ni{k) = -4 (nl{k) -n l d {k) 

4(m + uj + k z f - k\ 
9 (m + uj + k z ) 2 + kf 



n{k), (9) 



one can check that the total angular momentum sum rule is automatically 
fulfilled at the hadronic scale: 
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(10) 



Another interesting relationship connects L z to the longitudinal (gi) and 
transversity (hi) parton distributions, namely [0: 

gt(x,v%) + L a z (x,tf) = hZ(x,VLl), (11) 

and is naturally fulfilled in our approach. This relationship also holds for other 
relativistic models of the nucleon, such as the bag model. Let us stress that 
Eq. (PD is valid at the hadronic scale /Xq only and one should be careful when 
using it to extract information about L a z because it is broken by evolution, 
even at small values of Q 2 . This can be easily demonstrated by considering 
the singlet combination corresponding to Eq. ([11]), i.e. 

E{x, fil) + L z (x, f4) - H{x, nl) = (12) 

where S(x,/xg) = E a (^(x,/^) + gf(x,fx 2 )) and H(x,f4) = E a (/iJ(x,^g) - 
hi(x, /io))Q- m order to check the validity of Eq. (jTTf) at Q 2 > /Iq let us evolve 
(at LO) the first moments of the left-hand side of Eq. fll2|): 



(Y,(x, Q 2 ))i + [L z (x, Q 2 ))i - (H(x, Q 2 ))i = - r 50 / 81 )(E(x, fi)) x 

+ (b-^-b-^)(H(x,^h 

+ ^(i-^ 50781 ) (13) 

where b = ^^^pj- Clearly, the right-hand side of the equation above vanishes 
only if Q 2 = /1q. Furthermore, due to the form of the 6-dependent coefficients, 
it quickly deviates from at the initial stages of evolution pointing out the 



limits of the attempt, carried out in [|Tj], of extracting information on L z from 
Eq. flU]). 

Coming back to our evaluation of OAM, Eqs. (0) - @ show that the 
exact ratio between the amount of OAM and spin will depend on the specific 
form of n(k), or equivalently, on the spatial nucleon wave function. Let us 
note however that the momentum density averages many details of the spatial 
wave function and to this respect, the sensitivity of the final results to the fine 
details of the spatial wave function is reduced. 

-'^The minus sign in front of h\ comes from the properties of the operator that defines 
the transversity under charge conjugation, and therefore the analogous of Eq. (|ll]) for 
antiquarks should read g\{x, itg) + L a z {x,p^) — — hf (x, /j,q). Though our model does not 
contain antiquarks at the scale /Uq, this relationship can be easily checked in the bag model. 
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In the following we will discuss predictions obtained solving explicitly the 
mass equation: 



Mm = y K 2 + m2 + v ) ^ = E ^ ( 14 ) 

h 

with an hypercentral phenomenological potential: 

K=-~ + k,£ + A , (15) 

where £ is the hyperradius defined in the usual way and r, K\ and A are 



free parameters that are fixed by spectroscopy requirements fl6| , |23| , [24f . It is 
worthwhile mentioning that MR has no effects on the energy levels of the con- 
fining mass operator (|14]) - ([15]) explaining to some extent the success of non- 
relativistic (or relativized) approaches in reproducing the baryonic spectrum. 
On the other hand, another remarkable effect of the relativistic mass equation 
is the enhancement of the high momentum components in the nucleon wave 
function. Since the MR factor involves momentum dependent terms, the final 
results will be biased by the presence of these high-momentum components. 
In order to test the sensitivity to the details of the momentum density we will 
consider an additional scenario where the MR factors are combined with a 
wave function obtained from the non-relativistic Schrodinger reduction of the 
Eq. (|14"D with the same form of potential (|T5|). This new spatial wave function, 
hereafter indicated by ty', will contain far less high momentum components. In 
fact one of the risks of guessing the wave function instead of solving the mass 
equation (|14"D explicitly, is to underestimate the contribution coming from the 
high-momentum components of the correct solution, mostly carried over by the 
relativistic kinetic energy operator in the mass equation. Although the use of 
MR is not fully consistent when fy' is considered, since it was derived from a 
non-relativistic mass equation, we will discuss it as a 'pedagogical' example 
that represents an extreme scenario where high-momentum components have 
been strongly suppressed. The comparison of results obtained with \l/ and 
will serve to establish bounds on the effects of MR. 



3 Results and discussion 

The obtained OAM distribution at the hadronic scale /Zq, Eq. (^), for the 
wave function the solution of Eq. (|i"4"D, is shown in Fig. l.a. The outcome 
for the modified scenario (corresponding to the wave function \1>') is shown in 
Fig. l.b to appreciate the effect of the lack of high momentum components. 
Furthermore, the comparison with the bag model results [0 is also provided 
(Fig. l.c). It is clear that the LFCQM, regardless of the details of the spatial 
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wave function, provides OAM distributions which are comparable (even bigger 
by a factor 2) to the bag model. From the comparison between Figs, l.a and 
l.b one can see that the MR (and not the specific shape of the spatial wave 
function) is responsible for this sizeable OAM. In non-covariant quark models 
such as the Isgur-Karl model, where MR is omitted, the OAM distributions is 
almost flat pi. Even when considering a D-model P5[ where the probability 



of the D-wave component is raised up to a 20 %, the resulting OAM, though 
comparable in size to those obtained here, are peaked at lower x. Nonetheless 
the large deformation of the nucleon in the D-model should not be taken as 
realistic. 

In order to bring the OAM distributions to the high-energy experimental 
scale, we use the recently obtained evolution equations at LO [§, p6| . In the 
process the OAM distributions for the gluons will be generated. The initial 
scale /Xq is determined following the criteria exposed in |L6|], and at LO turns 
out to be /Zq = 0.079 GeV 2 . In fig. l.a and l.b we also present the evolved 
OAM distributions up to Q 2 = 10 GeV 2 (short-dashed line) and Q 2 = 1000 
GeV 2 (long-dashed line). 

By comparing again the LFCQM with the bag model (Fig. l.c) it is clear 
that a non- vanishing OAM persists in the large x region and this is a distinctive 
feature of relativistic treatments of the nucleon. Indeed, in I-K models, the 
OAM is entirely concentrated at low x. This may constitute a clear signature 
of relativity in the low-energy models of the nucleon if L z (x, Q 2 ) is measured. 

In our approach all the gluon OAM is generated through evolution. In 
Fig. 2 we present the resulting L g (x,Q 2 = 10 GeV 2 ) for \I>, and the bag 
model. There is an inverse correlation between the amount of high momentum 
components in the wave function and the value of L g at small x. The OAM 
gluon distribution for the bag model falls between those obtained with \1/ and 

Concerning the first moments of the distributions, our model gives a value 
for L q (fil) = J L z (x,fil) dx that ranges from 0.272 to 0.126 (\I/ and model 
respectively). It should be stressed that the corresponding values for AE (0.456 
and 0.748 respectively) are per se a clearcut signature in favor of light-front 
quark models, when compared to recent analysis of data (AS = 0.45 ± 0.09) 
II . Furthermore, these numbers are quite close to the angular momentum 
share-out given by the bag model. 

The first moments that make up the spin sum rule also evolve with Q 2 



according to 26 



Lq{Q 2 ) 



(6"i - l)^AE(^) 



+ b 8iL ? (/^ 



— (6 si 
50 v 



(16) 
(17) 
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J g (Q 2 ) = b-*J B Q4)--(b-*-l) (18) 

In Fig. 3 we show the evolution of these quantities with Q 2 for \1/ (Fig. 
3. a) and ty' (Fig. 3.b). It is worthwhile mentioning that, even if we do not 
have gluons at the hadronic scale, they quite rapidly develop a sizeable angular 
momentum content. Furthermore the gluon angular momentum evolves decou- 
pled from the quark sector and if we start with a vanishing J g (/J%) then J g {Q 2 ) 
is completely determined by the QCD anomalous dimensions. The values for 
J g (Q 2 ) in the region between 1 and 10 GeV 2 that we find (J g ~ 0.20 - 0.25) 
are compatible with those found by using QCD sum rules p7| (J g ~ 0.25) 



and in a recent lattice calculation |28[ (J g = 0.20 ± 0.07). They also agree 
with another model calculation based on the one-gluon exchange interaction 
between quarks |29| (J g ~ 0.24). The consideration of a non-vanishing J g ($) 
would not change much our results since it would also raise the scale /x 2 , (due 
to the fact that at that scale the gluons would carry some momentum) and 
hence h would be larger for a given Q 2 . 

Though the large error bars in the first direct measurement of the ratio 
As A? = o.41 ± 0.18(stat) ± 0.03(syst), and the values for Ag obtained 
in recent data analysis {§], Ag(Q 2 = 1 GeV 2 ) = 1.6 ± 0.9, do not allow to 
discriminate between models, our results fall within the range of the latter. 
As a matter of fact the rather moderate values for J g result from a strong 
cancellation between L q (Q 2 ) and Ag{Q 2 ) and, in particular we have Ag{Q 2 = 
1 GeV 2 ) = 1.36 and 2.22 for * and respectively. 

In a non-relativistic quark model one would expect L q {Q 2 ) ~ —J g {Q 2 ) ~ 
—0.25 in the range Q 2 ~ 1 — 10 GeV 2 since AS is a constant with Q 2 . When 
relativistic spin effects are taken into account, as Fig. 3 shows, one expects 
L q {Q 2 ) to be much smaller (L q (Q 2 ) ~ —0.12 at most) or close to zero. 



4 Concluding remarks 

In summary we have shown that covariant light-front based quark models 
give rise to non-trivial predictions for the OAM distributions at both low and 
high momentum scales. This departure from traditional treatments of the 
angular momentum structure of the nucleon is more manifest in the high-x 
region of the quark sector. We have seen that the performance of LFCQM 
is quite similar to other relativistic models of the nucleon such as the bag 
model. This comparison holds for a quite flexible choice of the mass operator. 
We have studied the predictions for other potentials that interpolate between 
the two somehow extreme situations presented here and conclusions are not 
changed. In fact, there is a clear correlation between the amount of high- 
momentum components in the momentum density n(k) and the size of the 
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OAM distribution. A more realistic interaction would give results closer to 
those of ^ than to those obtained with *&' because a relativistic treatment of 
the kinetic energy operator inevitably emphasizes the high-momentum tail. 

One should keep in mind however that the origin of the relativistic aspects 
is not the same in the bag and in the LFCQM presented here. While in the 
former the non-vanishing OAM comes from the small Dirac components, in the 
latter these ones are absent and relativity enters through the momentum de- 
pendence of the Pauli spinors. Certainly other more sophisticated spin-flavour 
basis can be constructed, such as the Dirac- Melosh bases |3j] , where covariance 
is manifest. Despite the fact that the used basis contains only kinematic and 
not dynamical (higher Fock states) effects, it represents a minimal framework 
that combines in an elegant way simplicity and a proper treatment of boost. 
Results obtained with this basis and with the bag model are of similar quality 
pointing out that it allows an easy implementation of relativistic effects in the 
spin structure of the nucleon. 

We gratefully acknowledge Vicente Vento for useful comments and a careful 
reading of the manuscript. 
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Figure captions 

Figure 1 Quark orbital angular momentum distributions calculated in light- 
front dynamics with the wave function ^ (a), with the modified wave 
function (see text) (b) and in the bag model (c). Solid lines correspond 
to the initial hadronic scale /z 2 ,, short-dashed lines to Q 2 = 10 GeV 2 and 
long-dashed ones to Q 2 = 1000 GeV 2 . 

Figure 2. Gluon orbital angular momentum distributions calculated at Q 2 = 
10 GeV 2 with the wave function \1/ (solid line), (long-dashed line) and 
the bag model (short-dashed line). 

Figure 3. The contributions to the proton spin sum rule according to the 
model with \1/ (a) and (b). The dashed curve shows |AS(Q 2 ), the 
long-dashed one is L q {Q 2 ) and the dot-dashed curve represents J g (Q 2 ). 
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